EPJ manuscript No. 

(will be inserted by the editor) 



Hypernuclei and in-medium chiral dynamics 

- - - Paolo Finclli a 

t> , 

' Physics Department, University of Bologna, 

O Via Irnerio 46, 40126 Bologna (Italy) 

<N 

CD ■ Abstract. A recently introduced relativistic nuclear energy density functional, 

constrained by features of low-energy QCD, is extended to describe the structure 
of hypernuclei. The density-dependent mean field and the spin-orbit potential of 
a /1-hyperon in a nucleus, are consistently calculated using the SU(3) extension 
of in-medium chiral perturbation theory. The leading long-range AN interaction 
arises from kaon-exchange and 27r-exchange with a X'-hyperon in the intermedi- 
ate state. Scalar and vector mean fields, originating from in-medium changes of 
the quark condensates, produce a sizeable short-range spin-orbit interaction. The 

■ model, when applied to oxygen as a test case, provides a natural explanation for 
^ ' the smallness of the effective A spin-orbit potential: an almost complete cancella- 
£h . tion between the background contributions (scalar and vector) and the long-range 

terms generated by two-pion exchange. 

7-H ■ 

qq ■ 1 Hypernuclear phenomenology 

| A hypernucleus is a nucleus in which one or more nuclcons have been replaced by strange 

■ baryons |l|2j . In particular, /1-hypernuclei are quantum systems composed of a single yl-hyperon 
plus a core of nuclcons. The standard notation, used throughout this work, is 



ON 

o 

: a 
o ■ A 



- 1—1 



total number of baryons (nucleons + hyperon) 
Z where ^ Z : total charge (not necessarily the number of protons) 
hyperon (in this case, but in general, also S, S, . . .) 



For instance, \ 3 C = 6p + 6n + A. A simple way to describe hypernuclear properties is to employ 
a phcnomenological potential [3] 



C/» = [/ c »+^(r), (1) 



where 



1 df(r) 



U?(r) = -V c A f(r) and U£ = V£( ) ^fl, (2) 

ym^c J r ar 

are the central and spin-orbit terms, respectively. f(r) is a radial form factor that can be 
determined from the corresponding nuclear density distribution f(r) — p(r)/p(0), or chosen in 
a Woods-Saxon form f(r) = 1/{1 +exp[(r — R)/a]}. In Fig.[T]we display the A binding energies 
Ba for a set of hypernuclei, as functions of A -3 / 2 . The dashed curves are extrapolations to 
the nuclear-matter limit (A — oo). When compared to the binding of a nucleoli, one notes a 
considerable reduction of Ba- The strength of the central potential V A is estimated 3j3 
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1 This analysis is based on data for light and medium hypernuclei (from \ 2 C to A°Ca), but the 
inclusion of heavier hypernuclei does not change the estimate of Eq. (|3]) significantly (28 MeV for a fit 
with a Woods-Saxon potential [3]). 
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Fig. 1. Data on binding energies Ba of A single-particle states (in various orbits) as functions of 
A~ 2 ^ 3 , where A is the mass number of the nuclear core (see also Ref. [4]). The experimental values 
have been taken from Ref. [5] (Tables 11 and 13). The dashed curves extrapolate the values of Ba to 
the nuclear-matter limit (A~ 2 ^ 3 — > 0). 



V C A ~ 32 ± 2 MeV , (3) 

i.e. about 2/3 of the depth of the nucleon potential. On the other hand, the A- nucleus spin-orbit 
interaction is extremely weak [3]: 

V£ ~ 4 ± 2 MeV . (4) 

This peculiar property has recently been confirmed by the E929 BNL experiment [6], which 
observed the spin-orbit splitting between {P\/2)a and (jp$/2)a orbits in JfC, by separately de- 
tecting the A inter-shell transitions (pi/2)^ — * (51/2)71 and (^3/2)71 — ► ( s i/2)a around 11 MeV 
(see Fig.HJ. The observed energy spacing is E(l/2~) - E(3/2 ) = 152±54±36 keV. Measure- 
ments of such high precision are feasible because the spreading widths of the A hypernuclear 
states are extremely small [5]. 

Considering that the corresponding spin-orbit strength for the nucleon is an order of mag- 
nitude larger (V^ ~ 20 MeV), it is extremely important to understand the microscopic mech- 
anism at the basis of this unnatural property. 



2 In-medium chiral dynamics applied to hypernuclei 

A novel approach to the nuclear many-body problem, based on in-medium chiral effective field 
theory, has recently been succesfully applied to nuclear matter and properties of finite nu- 
clei [7|8|9|10j . In particular, it has been demonstrated that iterated one-pion exchange and 
irreducible 27r exchange processes, with inclusion of Pauli blocking effects and Z\-isobar excita- 
tions in intermediate states, generate the correct nuclear binding |7l8j . The nuclear spin-orbit 
potential, on the other hand, is produced by the coherent action of scalar and vector mean 
fields, representing the in-medium changes of the quark condensates |9I10I11I12I13I15] . The im- 
portance of correlated 2ir exchange for hypernuclei was already pointed out in Ref. [16], and 
recently the in-medium chiral approach has been extended to include the strangeness degree 
of freedom [17j . The long-range AN interaction arising from kaon and 27r exchange, with a 
Z'-hyperon and medium insertions in intermediate states, has been explicitly calculated in a 
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Fig. 2. The hypernuclear 7 transitions observed in \ 3 C [516] . The levels are denoted by their excitation 
energy (in MeV), and the spin and parity. The observed spin-orbit splitting is small (154 ±54 ±36 keV), 
i.e. 20 ~ 30 times smaller than that of the nucleon p-levels. The ordering of the A p-levels appears to 
be the same as for the nucleon levels. (From Fig. 47 in Ref. [5]). 



controlled expansion in powers of the Fermi momentum kf. In Ref. [18) this approach has been 
applied to finite hypernuclei. 



2.1 The model 

In order to describe hypernuclei we extend the relativistic nuclear energy density functional 
(see Sec. 2.2 of Ref. [TO]), by adding the hyperon contribution: 

E [p} = E^[p]+E A [p}, (5) 

where Eq[p] describes the core of protons and neutrons (cf. Eq. (12) in Ref. pill), and E A [p] is 
the leading-order term representing the single /l-hyperon, decomposed in free and interaction 

E A [p}=E£ ee [p]+E? nt [p], (6) 



J d 3 r(M^AHl ■ V ± M a ]1>a\<M (7) 
J d\{(UG A s { P ) (U) (Ma) \M+ 

(<I> \Gv(p) (M) §A-f1>A) \<M) ■ (8) 

Here \<j)o) denotes the (hypernuclear) ground state. E^. ee is the contribution to the energy from 
the free relativistic hyperon including its rest mass Ma- The interaction term E A lt includes 
density-dependent hyperon-nuclcon vector (Gy) and scalar (Gg) couplings. They include mean- 
field contributions from in- medium changes of the quark condensates (identified with superscript 
(0)), and from in- medium pionic fluctuations governed by two-pion exchange processes (with 
superscript (7r)): 

G A (p) = Gf {a) +Gf M {p) with i = S,V. (9) 



with 



free 



El 
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Minimization of the ground-state energy leads to coupled relativistic Kohn-Sham equations for 
the core nucleons and the single /1-hyperon. Using the notation of Ref . [10] , they read: 

[-ry • V + M N + 70 (XV + Xfl + t 3 S tv ) + £ s + t 3 S T s} ipk = tki>k (10) 

[-*7 • V + M A + tqE§ + Si] i> A = tAipA , (11) 

where ipk and ipA denote the single-particle wave functions of the nucleon and the A, respec- 
tively. The single-particle Dirac equations are solved self-consistently in the "no-sea" approx- 
imation [19j . It is important to note that the rearrangement self-energy Er |20j is confined 
to the nucleon sector because all the density dependent couplings are polynomials in kf (and 
consequently in fractional powers of the baryon density through the relation p = 2 fc^/(37r 2 )), 
and there is no hyperon Fermi sea. The A self-energies read 

2# = C#(p)p, ££ = Gj(p) PS} (12) 

expressed in terms of the nuclear baryon and scalar densities, p and ps- 

In the following sections we analyze separately the different contributions to the density- 
dependent A- nuclear couplings Gf(p), arising from the kaon- and two-pion-exchange induced 
/1-nucleus potential, the condensate background mean-fields, and the pionic yl-nucleus spin-orbit 
interaction. 



2.2 Kaon- and two-pion exchange induced mean-field 

The density-dependent self-energy of a zero-momentum yl-hyperon in isospin-symmetric nuclear 
matter has been calculated in Ref. [17] at two-loop order in the energy density. This calcula- 
tion systematically includes kaon-exchange Fock terms (first diagram in Fig. [3]), and two-pion 
exchange with a X'-hyperon and including Pauli blocking effects in intermediate states. 



Fig. 3. One-kaon exchange Fock diagram 
and two-pion exchange Hartree diagrams 
with a X'-hyperon in intermediate states. 
The horizontal double-lines represent the 
filled Fermi sea of nucleons in the in- 
medium propagator: (~/-p— Mjv)[i(p 2 — M%+ 
ie)" 1 - 2tvS{ P 2 - M 2 N )e{ Po )6{k } - |p |)] [17]. 
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The self-energy is translated into a mean field A nuclear potential UA(kf). A cutoff scale A ~ 700 
MeV (or cquivalently, a contact term) represents short-distance (high momentum) dynamics 
not resolved at scales characteristic for the nucleon Fermi momentum. The value of A is adjusted 
to reproduce the empirical depth of the A nuclear central potential. 

Following the procedure outlined in Appendix A of Ref. [10j . we determine the equivalent 

density dependent A point coupling vertices G^\p) and Gy^\p). For the nucleon sector of 
the energy density functional the parameter set FKVW [10] is used. In Fig. [4] (case a) the A 
single-particle energy levels are plotted for the closed core of 16 nucleons (8n + 8p) plus a single 
hyperon (\ 7 0). At this stage of the calculation the p-shell spin-orbit partners are practically 
degenerate, and the energies of the doublets are, by construction, close to their empirical values. 
Even the calculated energy of the s-state is realistic, although slightly too deep in comparison 
with data (4 = -12.42 ± 0.05 MeV for \ 6 [5]). 

Up to this point the in-medium chiral SU(3) dynamics (including K- and 27r-exchange) 
provides the necessary binding of the system, but no spin-orbit force. As it has already been 
shown in Ref. |10] , the inclusion of derivative couplings does not remove the degeneracy of the 
spin-orbit partner states. 
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Fig. 4. yl single-particle energy levels in ]\0: (a) the contribution of the single-particle potential with 
density-dependent coupling strengths determined in Ref. [17] by including chiral K- and 27r-exchange 
(see Sec. I2.2[l ; (b) the spin-orbit effect from in-medium quark condensates has been included, with 
X — 2/3 corresponding to the simple quark-model prediction [21] (see Sec. 12.3ft : (c) the result of 
the additional compensating effect from the second-order AN tensor force with intermediate £ (see 
Sec. [22}). 



2.3 Background scalar and vector mean-fields 



In contrast to the mean-field induced by kaon- and two-pion exchange, the condensate back- 
ground A self-energies Sy^ and ZJg^ produce a sizeable spin-orbit potential, analogous to 
the nucleon case |10j . Finite-density QCD sum rules predict moderate Lorentz scalar and vector 
self-energies for the A hyperon. Under some reasonable assumptions about the density depen- 
dence of certain four-quark condensates [22pl . one expects a reduction of the corresponding 
coupling^ 



G 



A(0) 

s,v 



(13) 



by a factor \i where G^y and G^ are the vector and the scalar couplings to nucleons arising 



'v w s 

from in-medium changes of the quark condensates, (qq) and (q^q). The values of G\ 
have been determined by fitting to ground-state properties of finite nuclei [10| , and found to be 
in good agreement with leading-order QCD sum rules estimates 



and Gf 



2 While is rather insensitive to the details of the calculation, Eg y "' can only attain realistic 
values if the four-quark condensate (qq) 2 PN depends weakly on the nucleon density, and the four-quark 
condensate {qq) PN {ss) PN has a strong density-dependence [22] . 

3 For simplicity the same reduction factor \ is used both for the scalar and vector self-energies. In 
general one could use two parameters (xs and xv), but the difference can be easily absorbed by the 
cut-off A. 
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In Fig. 2] (case b) we plot the A single-particle energy levels calculated with the inclusion 
of these scalar and vector mean-fields, using the quark model prediction for the reduction 
parameter % = 2/3 |21j . In this case the p-shell spin-orbit partners are no longer degenerate. 
The calculated spin-orbit splitting is of the order of ~ 2 McV. However, the choice \ = 2/3 
is a rather simplistic estimate. A detailed QCD sum rule analysis suggests a reduction to 
X ~ 0.4 — 0.5 jl 1117122 ] , and to even smaller values if corrections from in-medium condensates 
of higher dimensions are taken into account. We note that the A- nuclear spin-orbit force is 
evidently still far too strong at this level, just as in the phcnomcnological rclativistic "sigma- 
omega" mean-field models. 



2.4 /1-nuclear spin-orbit interaction from chiral SU(3) two-pion exchange 

The /1-nucleus spin-orbit interaction generated by the in-medium two-pion exchange AN in- 
teraction has been evaluated in Ref. [IT]. In the spin-dependent part of the self-energy of a A 
hypcron scattering in slightly inhomogeneous nuclear matter from initial momentum p — q/2 
to final momentum p + q/2, one identifies the spin-orbit term S^(kf) = ^U£(kf) c-(qxp). 
It depends only on known SU (3) axial- vector coupling constants and on the mass difference 
between the A and S. The relevant momentum space loop integral is finite, and hence model 
independent in the sense that no regularizing cutoff is required. The result, 

Ufiikf*) -15 MeV fm 2 at kf ] ~ 1.36 fm" 1 , 

has a sign opposite to the standard nuclear spin-orbit interaction^. This term evidently tends 
to largely cancel the spin-orbit potential generated by the scalar and vector background mean- 
fields. It is important to note that such a "wrong-sign" spin-orbit interaction (generated by the 
second-order tensor force from iterated pion exchange) exists also for nucleons [2J]. However, 
this effect is compensated to a large extent by the three-body spin-orbit force involving virtual 
/_A(1232)-isobar excitations [25] (see Fig. [5]), so that the spin-orbit interaction from the strong 
scalar and vector mean-fields prevails. For a /1-hypcron, on the other hand, the analogous 
three-body effect does not exist, and the cancelation is now between spin-orbit terms from the 
(weaker) background mean-fields and the in-medium second-order tensor force from iterated 
pion-exchange with intermediate S. The small U — A mass splitting, Ms — Ma = 77.5 McV, 
plays a prominent role in this mechanism. 



Fig. 5. Three-body diagram of two-pion exchange 
with virtual Z\(1232)-isobar excitation (left). For 
a nucleon it generates a sizeable three-body spin- 
orbit force of the "right sign". The horizontal 
double-line denotes the filled Fermi sea of nucle- 
ons. The analogous diagram does not exist for a 
/1-hyperon (right), simply because replacing the 
external nucleon by a /1-hyperon introduces as the 
intermediate state on the open baryon line a E- 
hyperon, for which there is no filled Fermi sea. 



N 



D 



In order to estimate the impact of this genuine "wrong-sign" yl-nuclear spin-orbit term, we 
introduce 

2r dr 



Recall that nuclear Skyrme phenomenology gives Uu{kf) = 3WVo/2 ^ 30 MeV fm 2 for the 
strength of the nucleon spin-orbit potential [23] ■ 
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with the form-factor determined by the normalized nuclear density profile /(r) = p(r)/ p(r = 0). 
The corrections to the A single particle energies €a are then evaluated in first-order perturbation 
theory: 

e' A = e A + {<t>\AHi\<t>) , (15) 

where \<f>) denotes the self-consistent solution of the system of Dirac single-baryon equations 
(flC)]) and (fTTj) . In Fig. [4] (case c) one observes that the resulting p-shcll single-particle energy 
levels, corrected according to Eq. (fT5|) . are close to being degenerate. The spin-orbit splitting is 
now strongly reduced, but still rather too large in comparison with empirical estimates. This 
could be a consequence of the possibly too large quark- model reduction factor x = 2/3. We note 
that within the range of values of x compatible with QCD sum rules estimates, the empirical, 
almost vanishing spin-orbit splitting for the single- A states can indeed be obtained [TS] , 

In Fig. [HI we plot the A spin-orbit spacing 8 a = eA^Pi^) ~ e ^(^P3^ as a function of the 
ratio x between the background mean-fields for the yl-hyperon and for the nucleon. The cir- 
cles denote the spin-orbit splittings produced by the scalar and vector background fields alone. 
Even for unnaturally small values of x, the splitting remains systematically too large in com- 
parison with empirical estimates. Introducing the model-independent spin-orbit contribution 
from second-order pion exchange fEq. (fl~5|) ). these values are systematically reduced by about 
1.3 MeV (triangles). For x in the range 0.4 — 0.5 determined by the QCD sum rule analysis of 
Refs. |11I22) . the small spin-orbit splitting is now reproduced in agreement with recent empirical 
values [5]. 




Fig. 6. Evolution of the spin-orbit splitting 
5a = e A {lp~J 2 ) - e A (lp~^ 2 ) in \ 7 0, as a func- 
tion of the ratio x between the background self- 
energies of the A and nucleon. The dashed line 
at x — 2/3 denotes the simple quark-model 
value. Also indicated is the x-interval allowed by 
the QCD sum rule analysis of Refs. [11117122] . 
Calculations with (without) the chiral SU(3) 
spin-orbit correction (Eq. ([15}) are denoted by 
triangles (circles). The shaded area represents 
an estimate of 8 A (-0.8 MeV < 8 A < 0.2 
MeV) based on the measured energy difference 
AE{2+ - Of) in \ g O [5]. 



3 Conclusions 

The compensating mechanism for the spin-orbit interaction of the yl-hyperon in nuclear mat- 
ter, suggested in Ref. |17j . successfully explains the very small spin-orbit splittings in finite 
A hypernuclei. We emphasize that this mechanism, driven by the second-order pion-exchange 
tensor force between A and nucleon, with intermediate Z'-states, is model-independent in the 
sense that it relies only on SU(3) chiral dynamics with empirically well determined constants. 
This intermediate-range effect (independent of any regularization procedure) counteracts short- 
distance spin-orbit forces. In ordinary nuclei, i.e. without hyperons, the corresponding effect 
is neutralized by three-body spin-orbit terms (induced by two-pion exchange with virtual A- 
isobar excitations). These terms are absent in hypernuclei. The theoretical framework based on 
in- medium chiral dynamics, described in this work, will be systematically applied to heavier 
hypernuclei. 
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